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ABSTRACT 

We report deep K-band (18-27 GHz) observations with the 100-m Green 
Bank Telescope of HCN(1 — 0) line emission towards the two submillimeter- 
selected galaxies (SMGs) SMM J02399— 0136 (z = 2.81) and SMM J16359+6612 
(z = 2.52). For both sources we have obtained spectra with channel-to-channel 
rms noise of o < 0.5 mJy, resulting in velocity-integrated line fluxes better than 
< O.lJykms -1 , although we do not detect either source. Such sensitive obser- 
vations - aided by gravitational lensing of the sources - permit us to put upper 
limits of ^hcn ~ 2x 10 10 Kkms -1 pc 2 on the intrinsic HCN(1 — 0) line luminosities 
of the two SMGs. The far-infrared (FIR) luminosities for all three SMGs with 
sensitive HCN(1— 0) observations to date are found to be consistent with the tight 
FIR-HCN luminosity correlation observed in Galactic molecular clouds, quiescent 
spirals and (ultra) luminous infrared galaxies in the local Universe. Thus, the 
observed HCN luminosities remain in accordance with what is expected from the 
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universal star formation efficiency per dense molecular gas mass implied by the 
aforementioned correlation, and more sensitive observations with today's large 
aperture radio telescopes hold the promise of detecting HCN(1 — 0) emission in 
similar objects in the distant Universe. 

Subject headings: cosmology: observations — galaxies: high-redshift - galaxies: 
ISM - galaxies: starburst 

1. Introduction 

The search for molecular lines at high redshifts (z > 1) offers one of the richest and 
most exciting avenues to follow in the study of galaxy formation and evolution (Solomon 
& Vanden Bout 2005). The detections of rotational transitions of CO in extremely distant 
quasars (Carilli et al. 2002; Walter et al. 2003), high-z radio galaxies (Papadopoulos et al. 
2000; De Breuck et al. 2005) and submillimeter (submm) selected galaxies (Frayer et al. f 998, 
f999; Neri et al. 2003), have established that the most luminous and extreme objects in the 
early Universe harbour vast amounts of molecular gas (~ f 10 ~ n M Q - Neri et al. 2003; Greve 
et al. 2005) and have large dynamical masses (~ f 11 M© - Genzel et al. 2003; Tacconi et al. 
2006). Modelling of the bulk physical conditions of the molecular gas in high-z QSOs have 
even been possible in cases where several CO lines have been detected (WeiB et al. 2005). 
The success of CO observations is primarily due to its much larger abundance compared 
to other species as well as its easily excitable low- J (CO J + 1 — > J, J < 3) rotational 
lines, which makes it an excellent tracer of the total amount of diffuse (n(H 2 ) ~ I0 2 ~ 3 cm~ 3 ) 
molecular gas. However, the diffuse gas only serves as a reservoir of potential fuel available 
for star formation, and is not actively involved in forming stars. 

In our own Galaxy, the sites of active star formation coincide with the dense cores of 
giant molecular clouds (GMCs). This dense gas phase (n(H 2 ) > I0 4 cm~ 3 ) is best traced 
using molecules with high critical densities such as HCN. This basic picture also appears to 
hold in nearby galaxies, as suggested by HCN(f — 0) surveys of local (z < 0.1) Luminous 
Infrared Galaxies (LIRGs - L m ~ I0 11 L Q ) and Ultra Luminous Infrared Galaxies (ULIRGs 
- L m ~ I0 12 L Q ), which found a remarkably tight correlation between IR luminosity and 
HCN line luminosity (Solomon et al. 1992; Gao & Solomon 2004a,b). Recently, this relation 
was shown to extend all the way down to single GMCs, thus extending over 7-8 orders of 
magnitude in IR luminosity (Wu et al. 2005). This was interpreted as suggesting that the 
true star formation efficiency (i.e. the star formation rate per unit mass of dense gas) is the 
same in these widely different systems. 



- 3- 



Attempts have been made at extending the correlation between IR and HCN luminosity 
to primordial galaxies at high redshifts (Barvainis et al. 1997; Solomon et al. 2003; Isaak 
et al. 2004; Carilli et al. 2005), resulting in four objects detected in HCN out of 8 objects 
observed. So far the results have been consistent with the local correlation, however, the 
small number of sources precludes a statistically robust conclusion as to whether the relation 
also holds at high redshift. Furthermore, the majority of high- z objects observed have been 
extremely luminous QSOs and HzRGs dominated by Active Galactic Nuclei (AGN). 

The detection of HCN(1 - 0) in the z = 2.286 QSO F 10214+4724 (Vanden Bout, 
Solomon & Maddalena 2004), which marked the first high- ,2 detection of HCN using the 
Green Bank Telescope (GBT), motivated us to test the capability of the GBT to probe 
the dense gas in a population of high-z galaxies less extreme, i.e. less AGN-dominated than 
F 10214+4724, but more representative of starburst galaxies in the early Universe. The goal 
of our study was to observe two of the most prominent examples of submillimeter-selected 
galaxies (see Blain et al. 2002 for a review) - a population likely to constitute the progenitors 
of today's massive spheroids, and thus fundamental to our understanding of galaxy formation 
and evolution. 

Throughout this paper we adopt a flat cosmology, with Vl m = 0.27, JIa = 0.73 and 
H = 71kms- 1 Mpc- 1 (Spergel et al. 2003). 

2. Sources 

The two submillimeter galaxies targeted with the GBT were SMM J02399— 0136 (z = 
2.808) and SMM J16359+6612 (z = 2.517) 1 . These two sources were chosen by virtue of 1) 
their gravitational lensing amplifications making them amongst the brightest submm sources 
known (Smail et al. 1997, 2002; Kneib et al. 2004), and 2) their strong CO line emission 
(Frayer et al. 1998; Genzel et al. 2003; Sheth et al. 2004; Kneib et al. 2005). Moreover, while 
J02399 has a moderate amplification factor (fi ~ 2.5) and thus is intrinsically very luminous, 
the extremely large amplification factor of J16359 (/i ~ 22) means that this source has 
about an order of magnitude smaller intrinsic luminosity. Thus, our observations could 
potentially explore differences in the dense gas mass fraction between massive and not so 
massive starburst galaxies at high redshifts. 



hereafter we shall refer to these three sources as J02399 and J16359, respectively 
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3. Observations and Data Reduction 

The observations were carried out between March 21st and April 5th 2005 at the 100-m 
Robert C. Byrd Green Bank Telescope in Green Bank, West Virginia 2 . Due to scheduling 
and observability constraints, J02399 received by far the most observing time: 9.6 hrs of 
effective integration in each polarisation, with J16359 getting only 3.3 hrs. 

At the redshifts of the two sources, the HCN(1 — 0) line (u ieat = 88.631 GHz) falls within 
the K-band (17.6-27.1 GHz) receiver on the GBT. The K-band receiver consists of two pairs 
of dual circular polarisation beams, covering the frequency ranges 17.6-22.4 GHz and 22.0- 
27.1 GHz, respectively. In our case, the HCN(1 — 0) line from both sources could be observed 
using only one pair of beams (22.0-27.1 GHz). At ~ 25 GHz the GBT has a beam size of 
fwhm = 30" - much larger than the source size (< 2"). 

The observations were made in dual-beam mode, in which one feed sees the source (ON) 
and the other sees the blank sky (OFF). After one such scan (which in our case lasted 2 
minutes), the telescope nods so that the first feed now is OFF and the second feed is ON 
the source, thus ensuring that one beam is always integrating on-source (except during the 
nodding). As backend we used the GBT Spectrometer in its moderate-bandwidth, low- 
resolution submode, giving 200 MHz of total instantaneous bandwidth divided into 8192 
channels each of width 24.41 kHz. At 25 GHz this corresponds to a velocity-coverage of 
~ 2400 kms -1 and a resolution of ~ 0.3 km s -1 per channel - more than enough to encompass 
the entire line and resolve it spectrally. This is based on the reasonable assumption that the 
HCN line widths are similar to those of CO (Greve et al. 2005). 

Calibration was by noise-injection, in which a noise-diode is switched on and off at a 
rate of 1 Hz during the integration. Thus a typical data time stream can be visualised as 
follows: 



Beam 1 Beam 2 



ONi i 0FF 2 i 



Scan 1: [0N C(l/on , 0N mio// ] [0FF ca/on , 0FF caio// ] 

0FFi j2 0N 2 , 2 

/ * * / * * 

Scan 2: [0FF C(l/on , 0FF caio// ] [0N ca / on , 0N caio// ] 



2 The GBT is operated by the National Radio Astronomy Observatory. The National Radio Astron- 
omy Observatory is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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0Ni >3 0FF 2 , 3 
, * * / * * 

Scan 3: [0N ca / on , 0N caio //] [OFF ca / on , OFF ca i of f] 



The system temperature on the blank sky is then given by: 

Tsys{V,P) = T cal {u, P )— + , (1) 

Ufi* colon — Urr C aZo// * 

where T ca i[y^ p) is the calibration temperature of the noise diode and is a function of frequency 
(z/) and polarisation (p). Instead of T C(l /(z/, p) we adopted the mean T ca i{p) value across the 
entire bandpass. Throughout the observing run the system temperature (integrated over the 
entire bandwidth observed) varied between 28 K in optimal conditions to 55 K in the worst 
cases, depending on elevation and the water vapour content of the atmosphere. Each scan 
was inspected by eye in order to weed out bad or abnormal looking scans. 

The data were reduced using two independent methods. The first method used the 
GBTIDL-routine GETNOD, which takes the OFF signal from one scan and subtracts it from 
the ON signal in the neighbouring scan, and vice versa. In this scheme the spectra are 
calibrated onto the antenna temperature scale following: 

T antJ (v, P ) = r^(^p) ONl - 3 n r OFFl ' 3+1 , (2) 

where 0N i;j = (0N mion + 0N caio// )/2 and 0FF iJ+1 = (0FF ca/on + 0FF caio// )/2, and i, j represent 
the beam and scan number, respectively. Thus two neighbouring scans result in two inde- 
pendent spectra, T aniJ (z/, p) and T antJ+ i(z/, p). All the difference spectra are then averaged, 
weighted according to their system temperature, to produce the final difference spectrum for 
each polarisation. 

As an alternative reduction scheme, we used the method of Hainline et al. (2006) which 
involves constructing a baseline template for each ON signal from a best-fit linear combination 
of its neighbouring OFF signals. We can write this as 



T ant , j+1 (u,p) = T' sys (u,p) i ' J+ n ^ tempi ;' J+i; , (3) 



W lt i+1 - 0FF% $) 

where T' = a ■ T syS;i j + b ■ T syS;i j +2 is a combination of the system temperature of the scans 
that make up the baseline template: OFF*™^ = a-OFFij + b-0FF±^ +2 . The parameters a and 
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b are the best-fit parameters which minimises the quantity ((ONij+i — OFFi^^/OFF* 6 ^) 2 . 
Again, a final spectrum for each polarisation is constructed as a weighted average of all 
difference spectra produced in this manner. Vanden Bout et al. (2004) proposed another 
way of removing problematic baselines, which involves fitting a baseline template to the 
time-averaged spectrum containing the line. However, from a detailed comparison between 
the two methods, Hainline et al. (2006) demonstrated that their scheme is just as good at 
removing baseline-ripples at the scales one would expect to see a galaxy emission line. As a 
result we did not attempt to apply the method by Vanden Bout et al. (2004). 

The final spectra were corrected for atmospheric attenuation using the opacity values at 
23-26 GHz for the observing dates (R. Maddalena, private communication), and converted 
from antenna temperature to Janskys by applying an appropriately scaled version of the 
GBT gain-elevation curve at 43 GHz. 



4. Results & Discussion 

The final HCN(1 — 0) spectra of J02399 and J16359, reduced using the two indepen- 
dent methods described in the previous section, are shown in Figure 1. The spectra have 
been smoothed to SOkms" 1 bins. Both methods fail to get completely rid of residual base- 
line wiggles in the final spectra. However, the spectra reduced with GETNOD appear to 
be slightly less noisy than the spectra reduced with the second reduction method. The 
channel-to-channel rms noise of the spectra reduced using the GETNOD routine are 0.1 and 
0.3 mJy, respectively. This is somewhat higher than the theoretical noise estimates of 0.07 
and 0.12 mJy, calculated for a typical system temperature T sys = 38 K and integration times 
corresponding to the ones of J02399 and J 16359, respectively. Thus, it would seem that 
although the noise does integrate down as t~ 1//2 , the residual baseline wiggles prevents us 
from reaching the theoretical noise limit. 

No emission is detected significantly above the noise at Vlsr — Okms -1 (or any other 
Vlsr), which corresponds to the CO redshift and thus the expected position of the HCN(1 — 
0) line. However, the sensitivity of the observations allow us to put upper limits on the 
HCN(1 — 0) line luminosity of these two SMGs. In doing so we use the upper line flux limits 
derived from the spectra produced by GETNOD, as it results in the lowest noise spectra. 
The 3-er upper limits are calculated following Seaquist, Ivison & Hall (1995) 

S HCN(1 . 0) AV < ZatfvAv^) 1 ' 2 , (4) 

where a is the channel-to-channel rms noise, 8v the velocity resolution and At> fwhm the line 
width. The spectra were binned to a velocity resolution of 50kms _1 . We set the HCN(1 — 0) 



-7- 



line widths equal to that of the CO lines - see Greve et al. (2005). This is probably a 
conservative estimate since in local ULIRGs, the HCN line widths are rarely larger than 
those of CO (Solomon et al. 1992; Gao & Solomon 2004a). The resulting upper line flux 
limits are 0.08 and 0.13 Jykms -1 for J02399 and J16359, respectively. From these flux limits 
we derive upper limits on the apparent HCN(1 — 0) line luminosities of Lg CN(1 _ 0) < 5.0 x 10 10 
and < 6.6 x 10 10 K km s _1 pc 2 for J02399 and J16359, respectively. Correcting for gravitational 
lensing we find intrinsic line luminosities of L' HCN(10) < 2.0 x 10 10 and < 0.3 x 10 10 K km s _1 pc 2 . 
These upper limits on the HCN line luminosity are similar to those achieved towards high- ,2 
QSOs using the Very Large Array (e.g. Isaak et al. 2004; Carilli et al. 2005). Table 1 lists 
our findings along with all high-z HCN observations published in the literature at the time 
of writing. It is seen that of 10 sources observed to date only 4 have been detected - the 
remainder being upper limits. 

In order to determine the star formation efficiency per dense gas mass, and the dense 
gas fraction, as gauged by the L Fm /L' HCN and L' HCN /L' co ratio, respectively (Gao & Solomon 
2004a,b), we need accurate estimates of the FIR and CO(1-0) luminosities. 

The FIR luminosities of SMGs are generally difficult to estimate accurately, largely due 
to the poor sampling of their FIR/submm/radio spectral energy distributions. Furthermore, 
in our case the situation is complicated by the fact that most, if not all, of the sources in 
Table 1 contain AGN, which may be at least partly responsible for heating the dust and 
thus powering the FIR emission. Unfortunately, FIR/submm data from the literature do 
not allow for a detailed modelling of a hot (> 100 K) dust component (heated by the AGN), 
and we are therefore unable to correct for any AGN contamination. The one exception is the 
z = 3.9 BAL quasar, APM 08279+5255, where the separate AGN vs. starburst contributions 
to the FIR luminosity have been determined by Rowan-Robinson (2000) who find an appar- 
ent FIR luminosity of L FIR ~ 1.0 x 10 14 L Q for the starburst (corrected to the cosmology 
adopted here). In the case of SMGs, extremely deep X-ray observations strongly suggest 
that while every SMG probably harbours an AGN, it is the starburst which in almost all 
cases powers the bulk (70-90 percent) of the FIR luminosity (Alexander et al. 2005), and 
any AGN contamination would therefore not dramatically affect our conclusions. While 
optical spectroscopy of J16359 shows no evidence to suggest the presence of strong nuclear 
activity in this source (Kneib et al. 2004), this is not the case in J02399, which appears to 
be a type-2 QSO judging from its X-ray emission and optical spectrum properties (Ivison 
et al. 1998; Vernet & Cimatti 2001). Thus, J02399 could be a rare case where the AGN 
contribution is substantial. However, given our inability to correct its FIR luminosity for 
AGN contamination, we estimate the FIR luminosities of both sources simply by fitting a 
modified black body with a fixed (3 = 1.5 to their rest-frame SEDs and integrating it over 
the wavelength range 40 — 120/xm, see Table 1. This is consistent with the way in which the 
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FIR-luminosities of the objects observed by Carilli et al. (2005) were calculated. The FIR 
luminosities for the remaining objects in Table 1 were taken from Carilli et al. (2005) and 
converted to the cosmology adopted in this paper. 

The CO luminosities were taken from the compilation of high- ,2 CO detections in Greve 
et al. (2005). These are mostly high- J CO detections (J = 3 — 2 or 4 — 3), and in order to 
derive the CO(l — 0) line luminosity, we assume optically thick, thermalised CO line ratios, 
i.e. (3 - 2)/(l - 0) and (4 - 3)/(3 - 2) ~ 1, see Table 1. It should be noted, however, that 
studies on the ISM in local starburst nuclei yield CO (3 — 2)/(l — 0) ~ 0.64 (Devereux et 
al. 1994), while first results on ULIRGs yield a CO (6 - 5)/(4 - 3) ~ 0.6 ratio but with 
CO transitions of J = 4 — 3 and higher tracing a different H 2 gas phase than the lower 
three (Papadopoulos, Isaak & van der Werf 2006). These results along with known examples 
of very low global (high— J) /{km— J) CO ratios in high-z starbursts (e.g. Papadopoulos & 
Ivison 2002; Greve et al. 2003; Hainline et al. 2006) suggest that we may underestimate the 
SMG CO(l — 0) luminosities when we assume line ratios of unity, especially if observed CO 
J + 1 — > J, J + 1 > 4 lines are used. 

In Figure 2 a) - c) we plot L PIR against L' HCN , L' HCN against L' co , and L PIH against 
L' UCN /L' co , respectively, for our two sources along with all other existing high-z HCN obser- 
vations to date. For comparison we have also plotted the sample of local (U)LIRGs (open 
circles) observed by Gao & Solomon (2004a,b), which defines the local L IR -L^ CN relation: 
L m /L' UCN = 900 L & (Kkms' 1 pc 2 )" 1 . Since we are interested in the correlations with FIR- 
luminosity, we carefully fitted modified black-body SEDs to the Gao & Solomon sample and 
derived FIR-luminosities by integrating the SEDs from 40 — 120/im. To this end we used all 
available FIR/submm data in the literature for each source. The resulting local FIR-HCN 
correlation is given by logL FIR = (0.97 ± 0.07) logL' HCN + (2.9 ± 0.5), and is shown as the 
solid line in Figure 2 a). The slope of this correlation is consistent with unity, suggesting 
that the local FIR-HCN correlation is linear (Carilli et al. 2005). In the following we shall 
use this correlation, along with its 1-a limits allowed by the fitting errors (dotted lines in 
Figure 2 a)), to discuss whether the SMGs and other high-z sources follow the local FIR- 
HCN correlation. Secondly, motivated by what is seen in local (U)LIRGs, we shall impose a 
minimum dense gas fraction on the high- ,2 population and from that derive a lower limit on 
their HCN-luminosity. Finally, we assume the SMGs follow the local FIR-HCN correlation 
(within 1-cr) and explore what lower limits can be put on their dense gas fraction. 

With only four HCN detections to date and little more than a handful of upper limits, 
including the ones presented here, it is difficult to determine whether the local L PIH — L' HCN 
relation extends to higher redshifts and luminosities. It is interesting to note in Figure 
2 a), however, that while the three SMGs with HCN observations (our two sources plus 
SMM J14011+0252) all are consistent within 1-a of the local relation extrapolated to higher 
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luminosities, 5/7 of the QSOs observed in HCN fall above the 1-a envelope. This may sug- 
gest that QSOs have higher FIR luminosities for a fixed HCN luminosity than the SMGs. 
Although, this is not statistically significant, it is consistent with the expected picture in 
which the AGN in QSOs contribute a higher fraction to the FIR luminosity than in SMGs. 
Correcting for the AGN contribution in QSOs would lower them from where they are cur- 
rently plotted, bringing them more in line with the local L FIR — L' HCK correlation. 

Alternatively, the data could in principle also be interpreted as a steepening of the 
slope of the correlation at higher luminosities, meaning that the AGN contribution in both 
QSOs and SMGs is higher than in local (U)LIRGs - a scenario which would be in line with 
the suggested increase in AGN dominance with FIR luminosity in local ULIRGs (Veilleux 
et al. 1999; Tran et al. 2001). Although, X-ray constraints suggest that AGN contribute 
< 30 percent of to the FIR luminosity of SMGs, they do not allows us to completely rule 
out the above scenario since it is possible many SMGs harbour obscured, Compton-thick 
AGN (Alexander et al. 2005). However, the fact that the star formation rates for SMGs 
derived from reddening-corrected Ha luminosities agree with those derive from their FIR 
luminosities, suggest that AGN are unlikely to contribute signficantly to the FIR. 

We see from Figure 2 b) that all the high- z objects lie above the L' HCN — L' co correlation 
fitted to local galaxies with FIR luminosities less than 10 11 L Q and with a fixed slope of 
unity. This is consistent with the the steepening of the L' HCN — L' co correlation at higher FIR 
luminosities seen within the local (U)LIRGs sample itself as reported by Gao & Solomon 
(2004b). They argued that the dense molecular gas fraction, as gauged by the L' UCN /L' co 
ratio, is a powerful indicator of starburst dominated systems. Also, locally they found that 
all galaxies with L' UCK /L' co > 0.06 have L m > 10 11 L Q . In Figure 2 c) we have plotted 
our estimates of L FIR for the local (U)LIRG sample against their L'^^/ L' co ratios, and find 
that apart from two sources, the same criterion can be applied to the FIR luminosity, i.e. 
sources with L' HCN /L' co > 0.06 have L PIR > 10 11 L Q . It is also seen that all the high- z sources 
observed so far - all of which have been very luminous - obey the same criterium. Given their 
large FIR luminosities, it therefore seems reasonable to assume that the SMGs too have high 
dense gas fractions, i.e. L' HCN /L' co > 0.06. In that case, one can put a lower limit on their 
HCN luminosity and investigate where they would lie in the FIR-HCN diagram. The range 
of allowed HCN luminosities derived in this way for each SMG (and other high-z objects) 
is shown as horizontal lines in Figure 2a. It is seen that while J02399 and J16359, and the 
other high-z objects, could lie above the upper dotted line and still have a high dense gas 
fraction (L' HCN /L' co > 0.06), this is not the case for J14011. Thus, if our above assumption 
is correct, we can conclude that at least one of the SMGs (J14011), and possibly more, is 
consistent with the local FIR-HCN correlation, and that this SMG should be detectable in 
HCN if the sensitivity is improved by a factor x2. 
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APM 08279+5255 has the highest measured dense gas fraction, based on detections, in 
the high-z sample, even though it falls above the local L FIR — L' HCN relation in Figure 2 a). 
Naively, this would indicate that the ISM in this system is dominated by dense gas, and 
that star formation contributes substantially to its FIR-luminosity. However, caution should 
be taken since we have calculated the HCN/CO ratio using the nuclear CO line luminosity 
of APM 08279 (Papadopoulos et al. 2001). Adopting the global CO line luminosity, which 
includes emission extended on 10-kpc scales (Papadopoulos et al. 2001), i.e. scales comparable 
to those probed by the HCN observations (fwhm ~ 7" - Wagg et al. 2005), results in a 
HCN/CO ratio almost an order of magnitude lower. Similarly, we see that the upper limits 
on the dense molecular gas fraction in J02399 and J 16359 are rather high, yet both sources 
lie close to the L FIH — L' HCN correlation in Figure 2 a). What would the dense gas fraction 
of the SMGs be if they were required to lie between the upper 1-a line of the local FIR- 
HCN correlation and their observed HCN limits (for unchanged FIR luminosities)? The 
answer is given in Figure 2 c), where the permitted range in HCN/CO luminosity ratios are 
marked as horizontal lines. It is seen that J14011 has a dense gas fraction well within the 
range of most ULIRGs, and could have L' HCN /L' co < 0.06 and still be consistent with the 
FIR-HCN correlation. For the other two SMGs we see that even in the most conservative 
case, where the SMGs are just barely consistent with the 1-a line, they would still have 
some of the highest dense molecular gas fractions observed, comparable to the most extreme 
starburst dominated ULIRGs seen locally. Thus we conclude that if the SMGs do in fact 
follow the FIR-HCN correlation (or are within 1-a of it), then J02399 and J 16359 are likely 
to have higher dense gas fractions than J14011. Note, however, that if an AGN contributes 
significantly to the FIR luminosity, the above argument would yield a higher HCN luminosity, 
and thus higher HCN/ CO ratios than its true value, but we emphasize that the energy output 
at FIR/submm wavelengths from AGN in SMGs is found to be negligble compared to the 
starburst (Alexander et al. 2005). 



5. Summary 

With this study we have increased the number of submm-selected galaxies with HCN(1 — 
0) observations from 1 to 3, although the upper limit on the HCN luminosity of one of our 
sources (J16359) is not a very sensitive one. The upper limits on the HCN(1 — 0) line 
luminosity of these three SMGs are consistent with the L FIR — L' HCN relation observed in the 
local Universe (Gao & Solomon 2004a,b). In order to test whether a typical SMG with a 
FIR-luminosity of L FIR ~ 10 12 ' 5 L Q is inconsistent with the local relation, we need to be 
sensitive to HCN luminosities of L' HCN < 10 10 Kkms -1 pc 2 . Thus, we appear to be close (less 
than a factor of two) in sensitivity to that needed from Figure 2 in order to detect HCN 
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in SMGs. We conclude that until deeper observations come along there is no evidence to 
suggest that SMGs have higher star formation efficiencies per unit dense gas mass than local 
(U)LIRGs or indeed Galactic GMCs (Wu et al. 2005). 
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Fig. 1.— HCN(1 - 0) spectra of J02399 (left) and J16359 (right) produced using the 
GETNOD routine in GBTIDL (top panels) and our own reduction technique described in §3 
(bottom panels). The spectra have been plotted on the same flux- and velocity-scale, mak- 
ing it easier to compare between the different sources and reduction techniques. The LSR 
velocity-scale is relative to the CO redshift. The spectra have been binned to ~ 50kms _1 
and the channel-to-channel noise varies between O.f — 0.5 mJy 
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Fig. 2. — Correlations between different combinations of FIR, HCN and CO luminosity for 
local (U)LIRGs (open circles - Gao & Solomon 2004a,b) and the 10 high-z objects observed 
in HCN to date (solid circles and triangles). The latter have been corrected for gravitational 
lensing where appropriate, see Table 1. The three grey triangles denote upper limits of 
the two SMGs observed as part of this study and SMMJ14011+0210 which was oberved in 
HCN by Carilli et al. (2005) using the Very Large Array. The solid line in a) represents the 
local FIR-HCN correlation derived in this paper: logL FIR = 0.97 log L^ CN + 2.8, while the 
horizontal extensions to the upper limits mark the L' HCN -range allowed for each source given 
its observed upper HCN limit and assuming that L' HCN /L' co > 0.06. The solid line in b) is 
a fit (with a fixed slope of 1) to normal spiral galaxies in the local Universe. The dotted 
line in c) marks the L' UCN /L' co limit above which sources in the Gao & Solomon (2004a,b) 
sample are (U)LIRGs. The range in HCN/CO luminosity ratios for the SMGs, if they were 
to lie between the upper 1-a line of the local L FIR — L' HCN correlation, and their observed 
HCN limits are shown as grey horizontal lines. 



Table 1. Observations of HCN at high redshifts to date. The two SMGs studied here are in the upper part of the 
table. The luminosities have not been corrected for gravitational lensing, but the lensing amplifications (/i) are given 

in column 9. 



Source 


Type 


z 


HCN line 


ji 

HCN 


CO line 


r' 

^CO 






Ref. 










Xl0 10 Kkms-ipc 2 




xlO 10 Kkms- 1 pc 


2 Xl0 12 L Q 






SMM J02399-0136 


SMG 


2.808 


1 - 


< 4.9 


3-2 


12 + 2 


30 


2.5 


[1].[2],[3] 


SMM J16359+6612 


SMG 


2.517 


1-0 


< 6.7 


3-2 


8.4 + 0.4 


22 


22 


[1],[4],[5] 


SMM J14011+0252 


SMG 


2.565 


1-0 


< 0.8 


3-2 


9.0 + 0.7 


8 


5 


[6], [7], [8] 


SDSS J1148+5251 


QSO 


6.419 


1-0 


< 1.0 


3-2 


2.7 + 0.6 


28 


1 


[6], [9], [10] 


BR 1202—0725 


QSO 


4.693 


1-0 


< 4.4 


2-1 


9.2 ± 1.0 


62 


1 


[11], [12] 


APM 08279+5255 


QSO 


3.911 


5-4 


4.2 ± 0.5 


1-0 


14 + 3 


100 


3 a 


[13], [14], [15] 


MG 0751+2716 


QSO 


3.200 


1-0 


< 1-7 


4-3 


17 + 1 


20 


16.6 


[6], [16] 


QSO J1409+5628 


QSO 


2.583 


1-0 


0.7 + 0.2 


3-2 


8.3 ± 0.3 


34 


1 


[6],[17],[18] 


H1413+117 


QSO 


2.558 


1-0 


3.2 ± 0.6 


3-2 


43 + 7 


68 


11 


[19], [18], [19], [20], [21] 








4-3 


8.6 ± 1.6 










[22] 


IRASF10214+4724 


QSO 


2.286 


1-0 


2.2 ± 0.4 


3-2 


12 + 3 


48 


13 


[23], [24] 



References. — [1] This work; [2] Frayer et al. (1998); [3] Genzel et al. (2003); [4] Sheth ct al. (2004); [5] Kneib et al. (2004); [6] Carilli et al. (2005); [7] 
Fraycr ct al. (1999); [8] Dowries & Solomon (2003); [9] Bertoldi ct al. (2003); [10] Walter et al. (2003); [11] Isaak et al. (2004); [12] Carilli ct al. (2002); [13] 
Wagg et al. (2005); [14] Papadopoulos et al. (2001); [15] Lewis ct al. (2002); [16] Barvainis ct al. (2002); [17] Bcclcn et al. (2004); [18] Hainline et al. (2004); 
[19] Solomon et al. (2003); [20] WeiB et al. (2003); [21] Wilner et al. (1995); [22] Barvainis et al. (1997); [23] Vandcn Bout ct al. (2004); [24] Solomon ct al. 
(1992). 

a We adopt the lensing amplification derived from high-resolution CO maps by Lewis et al. (2002). Note, however, the amplification factor may be as high 
as fi ~ 7 - 20 (Downcs ct al. 1999) 



